NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2364 


ON TWO-DIMENSIONAL FLOW AFTER A CURVED STATIONARY 
SHOCK (WITH SPECIAL REFERENCE TO THE PROBLEM 
OF DETACHED SHOCKWAVES) 

By S. S. Shu 

Massachusetts Institute of Technology 



OObSbEt, 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECH URRARY KAFB, NM 




TECHNICAL NOTE 2364 

ON TWO-DIMENSIONAL FLOW AFTER A CURVED STATIONARY 
SHOCK (WITH SPECIAL REFERENCE TO THE PROBLEM 
OF DETACHED SHOCK WAVES) 

By S. S. Shu 

SUMMARY 


The problem of two-dimensional flow behind a curved stationary 
shock wave is considered analytically. The method assumes a given 
shock- wave shape, which automatically determines certain initial con- 
ditions on the flow variables; and the flow pattern, including any 
body shape, follows from the initial conditions. Approximate analytic 
expressions are found for the stream function in the subsonic region 
following the shock and, after the stream function is obtained, the 
flow density is determined by Bernoulli' s equation which connects the 
density with the derivatives of the stream function. The final solution 
can then be determined from the velocity field thus obtained. 


INTRODUCTION 


The problem of compressible flow after a curved shock wave has been 
investigated by various authors (references 1 to 4). In the analytical 
treatments, the main interest has so far been concerned with the local 
properties of the flow, such as the relations between the gradients of 
various physical and geometrical quantities along the shock wave and 
those along the body. 

The purpose of the present work is an attempt to treat analytically 
the two-dimensional problem after the shock in the large. The method 
assumes a given shock-wave shape, which automatically determines certain 
initial conditions on the flow variables. It is therefore a Cauchy 
problem and the flow pattern, including any body profile, follows from 
the assumed ini tial conditions. Approximate analytic expressions for 
the stream function are found for the compressible flow in the whole 
subsonic region. The stream function thus found satisfies the exact 
shock conditions but it satisfies the differential equation only approxi- 
mately. After the stream function is obtained, the flow density is 
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determined by the exact Bernoulli equation which connects the gas 
density with the derivatives of the stream function. Then the sonic 
line can be determined from the velocity field thus obtained. 

The author wishes to take this opportunity to t hank Professor 
C. C. Lin for his valuable criticisms and discussions. 

This work was conducted at the Massachusetts Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


ratio of specific heats of fluids 
free- stream pressure 
free-stream density 

free-stream velocity magnitude, assumed to be parallel 
to x-axis 

free-stream velocity of sound 
function determining shock shape (l) 
shock nose when F(o) =0 
stagnation density after shock wave 
stagnation velocity of sound after shock wave 
g^ - where ijr is stream function 

fl complex variable in holograph, plane 

((*x - if y)/{i p o^o 2 - S 2 + Po=o)j 

"f Q (r) parameteric representation of a three-dimensional curve 

with parameter r 

P,Q complex variables 


7 

Pi 

Pi 

U 

C 1 

F ( p l Uy ) 

(x o ,0) 

Po 

c o 
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x = — 5 F (pi u y) + x o 


equation for shock 


% 

t(x,y) 

tan t 
P 
P 

G(t) 

F(t) 

Z,m 

u,v 

q = \ju 2 + v 2 

* 

c local velocity of sound after shock wave 

maximum velocity magnitude of flow 

s,n tangential and normal directions at point on shock wave 

q^., q n tangential and normal components of velocity of flow 

immediately after shock wave 

Ki^Kg constants defined by equations (18) 

F(y) analytic function in real variable y 

0(x>y) auxiliary function defined by equations (22) 

h(ft) arbitrary analytic function of ft 

Re real part of a complex quantity 


free- stream Mach number (v 7 p i, 

stream function of flow after shock wave (ppUy on 
the shock) 

slope of shock wave at point on it 

pressure after shock wave 

density after shock wave 

entropy function 

function determining shock shape 

constants defined by equations (8) 

velocity components of flow after Bhock wave 




k 
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k = p^P/ PqCq 


z = x + yi 

f-^(t) analytic function for Weierstrass’ parametric represen- 

tation of a minimal surface 

t complex variable occurring in Weierstrass' formula for 

minimal surfaces 


Im 

^o 

I 

?(p),!(Q) 

k( r) 

*(y) 

tr = y + ico 
0) 


imaginary part of a complex quantity 
vector function defined "by equation (29) 

vector with, components x, y, and w 

complex-valued vector functions of variables P and Q, 
respectively; components denoted by subscripts 1, 2, 
and 3 

unit vector normal to tangential direction of ! 0 (r) 

vector function defined by equation (33) 
complex variable 

imaginary part of complex variable a 


FORMULATION OF PROBLEM 


The present report is concerned with a stationary shock wave pro- 
duced by a two-dimensional body placed in a ■uniform supersonic stream. 
The fluid is assumed to be an ideal gas and- the process on any material 
particle, to be adiabatic. Thus, the entropy is constant along a 
streamline while it may vary from one streamline to the next. 

Let the coordinate axes be so chosen that the incoming free stream 
is along the x-axis and that the shape of the shock wave is represented 
by 


- [f( p i Uy )/ p i D ] 


+ x_ 


X 


( 1 ) 
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where U is the free-stream velocity and p-]_ is the gas density in 
the free stream. As the shock wave must asymptotically become the 

Mach line at infinity, the function F 1 (p]Uy) — > tl^Mp^ - 1 as 
y — ^ ±co, where Mp is the free-stream Mach mmiber. 

Entropy of flow after shock waves .- The stream function ^(x^y) 
may he regarded as the amount of the flow across a curve joining the 
point (x,y) to the point (x o ,0). It must therefore he continuous across 

the shock wave because of the continuity of the mass. Thus, on the 
shock wave, 

iKx,y) = p-JJy (2) 


The slope of the shock wave is therefore 


tan r = ^ = l/F' (i|r) (3) 

dx 


Since the entropy is constant along a streamline, the pres sure- density 
relation for the flow after the shock wave can he written. in the form 


p/p 7 = gWv-l/p-J (4) 

where y is the ratio of the specific heats of the fluid, p-^ and p 
are, respectively, the pressures before and after the shock, and p-^ 
and p are corresponding densities of the fluid. This entropy function 
G(t 10 is now to he expressed in terms of the function F(tJ 0 which deter- 
mines the shock shape. From the shock conditions, it can he immediately 
shown that 

(5) 


(7 + l)Mp 2 sin 2 T 
2 + (7 r l)Mj_ 2 sin 2 T 



P/Pl - 


27M-J 2 sin 2 T + (l - 7) 


7+1 


(6) 
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Therefore 


G(t) = 


1 2rM x 2 + (1 - 7 ) [l + (F* )f] 


(7+l) r+1 1+(F«) 2 


-~|l + (F' )9 + 7 - 1 

M -, 2 


where F' = F'Cijr). Thus one has the formula 


(7) 


G' (i)/G(t) = 2F , (*)F"(*) 


my 


1+Z(F') 2 1 + m(F’ ) 2 1 + (F‘ ) 2 


(8a) 


where 


l = (1 - 7)/(27M 1 2 -7+l) 


(8b) 


m = 2 |[2 + (7 - 1 )M 1 ^] 


(8c) 


Fundamental differential equation for stream function .- By the 
definition of the stream function given above. 


PU 8y 




From the equation of vorticity (reference 3) 


5$ t.) + *r) - -®'(*)p r Pi/(r - 1 )p i 7 


and Bernoulli' s equation along a streamline 


q 1 - ___ c 1 - = c£ 7 P _ 

2 7 - 1 p 2 


i2 .2 

+ 

2 7-1 


(9) 


( 10 ) 


( 11 ) 
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it follows that 



2uv 


*: 


xy 



* 


yy 


= =i gKil 

27 G(t) 


2 2 

p 9 m 




.( 12 ) 


where 



+ v c 


(13) 


c 2 = 7P/p 


(14) 


and q^ is the maximum velocity of the flow. This equation is here- 
after referred to as the fundamental equation for compressible flow 
after the curved stationary shock wave given by equation (l). In it, 
the quantities u, v, c, and p are known functions of i|r, 
and iy determined by equations (9), (4), and (ll). 

It should be noted that equation (12) is equivalent to Crocco's 
equation (reference 5) for vortex motion. In fact, one is immediately 
derived from the other by the relation between the stream f unction 
defined above and that defined by Crocco. 


Shock conditions expressed in terms of derivatives of stream 
function . - Let s and n be, respectively, the tangential and the 
normal directions at a point on the shock wave. Then 


ty*. = ~ cos t + — sin r 


TX d E 


dn 


. . cM 

i|r v = ~ sin t - cos t 
J os on 


(15) 


From continuity of the flow across the shock it follows that 

= pq n = PpU Sin T 

cty TT 

~ = -pq^. = -pU cos r 


(16) 
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where q^ and q n are, respectively, the tangential and the normal 

components of the velocity of the flow immediately after the shock 
wave. Thus, by virtue of equations (15) and (6), the shock conditions 
in terms of the derivatives of the stream function are transformed into 
the following form: 


*x = Pl u 


F'(pi u y) . 1 + (f 1 ) 2 (pjUy) - 

1 + (P’) 2 (p 1 Uy) 1 + (p’) 2 ( Pl Uy) + LjJ: M-Z 


(f* ) 2 (p x uy) 1 + (f' ) 2 (p z Uy) - Mj 2 

Vy = Pl U - Pl U 

1 + (P> ) 2 (P 1 Uy) 1 + (p* ) 2 ( Pl Uy) + ini Mj 2 


(17) 


Mathematical' formulation of problem . - The problem of finding the 
two-dimensional-flow pattern after a stationary shock wave given by 
equation (l) is then to find the solution of the fundamental equa- 
tion (12) with the conditions (17) on the shock wave. It is therefore 
a Cauchy initial-value problem if the shock wave is given. 

However, it should be noted that the stream function Tjr(x,y) thus 
obtained is, in general, many-valued if it is extended continuously (in 
physical variables) into the whole region after the shock wave. In the 
actual case, it is expected that other shock waves are formed near the 
tail of the body. Therefore the solution, once obtained, is limited 
by such shocks. The flow farther downstream should be treated separately. 
In fact, the flow after the shock wave can be treated separately for the 
subsonic region and the supersonic region. 

Determination of physical variables .- After i}r(x,y) is found, the 
density distribution may be obtained from equations (4) and (ll). The 
velocity components may then be gotten from equation ( 9) . It should be 
remarked that this determination is two-valued: One value of the 

density corresponds to supersonic flow and the other, to subsonic flow. 

The location of the sonic line is therefore essential. When a solution 
is obtained for ijr (whether exact or approximate) the sonic line is 
immediately given by the condition that the two roots in Bernoulli' s 
equation for .the density p are equal. This gives 


( 


-LV " 1 
p iJ 


7 - 1/2 

7 + ly + 1 



= K 


1 


1 G(i|r) 
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Substituting back into Bernoulli' s equation, it follows that 

= K^rtj]^ 1 (18a) 

where 

*2 - <*»> 

and c-^ is the free- stream velocity of sound. 

In order that the sonic line thus determined should pass through 
the sonic point immediately after the shock wave, i x 2 +■ ijr 2 must 
satisfy the exact shock condition on the shock wave. 

Approximations . - In the following discussions, restrictions are put 
on the shock shape such that 

(1) The shock has its nose at (x Q ,o) where the shock is normal; 
that is, F'(0) = 0. 

(2) F(y) is a n alytic in the real variable y; moreover, F"(y) is 
positive. 

(3) The shock wave tends to .a Ifech line asymptotically. 

From the first two assumptions it follows that F'(y) is a monotonic 
function Of y and yF' (y) ^ 0. From the last assumption, it is seen 

readily that 0 ^F'(y) ^ - 1 j 1 ^ 2 . 

The effect of vorticity generated by the shock 'can be better under- 
stood by examining the term in the right-hand side of the fundamental 
equation (12) for i|r. From equations (8), the term on the right-hand 
side of equation (12) is a product of F"(i(r)p^ci M 2 and a nondimens ional 
quantity. By a straightforward calculation this nondimensional factor is 

(«i 2 -i ) 2 


7 


M ] 


2 + 2 1+ ) 2 
7-1 


f’U) 


1+ (F* ) 2 


1 - 


(F: )2* 


M-i - 1 


'1 + 
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A rough estimation shows that it is hounded by 


2 


1 F'(*) 1 1 

± (F’) 2 

2 (Ml 2 - l) 5/2 

(m x 2 - l) 1 / 2 1 + (F' ) 2 

M^ 2 - 1 

7 ^'\r )( Ml2 + r-i) 


It Is therefore clear that, for the free- stream Mach number close 

to 1, a perturbation process based on the parameter 



can he performed. The term in the right-hand side of equation (12) 
can he neglected for the first approximation. In fact, for any Mach 
number M^ in the general case,l the nondimens ional factor is an 

infinitesimal of first order for small values of ijf and of second 
order for large values of The assumption of neglecting the effect 

of vorticity is perhaps reasonable for a certain range of greater 

than unity. 

In the following discussion the effect of vorticity will he neg- 
lected in finding approximate solutions for the stream function. The 
subsonic region will now he treated analytically by finding approximate 
solutions for the stream function. The exact shock conditions will he 
satisfied hut the fundamental equation (12) will he satisfied only with 
the following approximations: 

(1) The effect of vorticity is neglected 

(2) Equation (12) is simplified by an approximate estimation of u/c 

and v/c as explicit functions and Tjy 

This estimation follows the same line as the one first suggested by 
Chaplygin (reference 6). The idea has been successfully applied by 
many authors (references 7 to 10) for analytic treatment of subsonic 
flows. 

^An upper bound for the nondimens ional factor in terms of M^_ only 
can be determined. The following is a trivial one: 
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It may further he suggested that the solution of the simplified 
equation for i}r satisfying the exact shock conditions may also serve 
as a first approximation for the supersonic region in the case of a 
thin body with a blunt nose. For such a case, the streamlines are 
likely to be nearly straight in the supersonic region. As the approxi- 
mate solution of ijr obtained from the above considerations represents 
the true geometrical pattern of streamlines on the shock wave, far away 0 
from the body, and near to it, it may be expected to be a reasonable 
approximation in the whole region. However, this can be applied only 
if the two-dimensional body is an analytic curve. In case the body has 
a discontinuity, the flow downstream should be treated separately from 
the Mach line passing through this point. 

After having obtained the first approximation (covering both the 
subsonic and the supersonic regions), an iteration procedure may then 
be carried out. This is done by evaluating all quantities in equa- 
tion (12) except j yy , ''J'xyj and tyyy from the approximate solution 

and then solving the linear equation for 

AN ANALYTIC APPROXIMATE SOLUTION IN THE LARGE 
FOR SUBSONIC RECTOR 


Solution of simplified differential equation .- The solution of 
equation (12) satisfying the initial conditions (17) on the shock will 
now be approximately obtained by finding the stream function ty(x,y) 
which satisfies the following simplified equation: 


I . ( *y f 

t X ty 

r, _/*fi 

\Po c o/ 

Vxx + p 0 *xy + 

Po c o 

0 

0 

0 

- 1 


*yy = 0 


(19) 


with the same initial conditions (17). Equation (19) can be written as 

= 0 ( 20 ) 


d 

Tjr ■~ 1 

y x 

+ af 

- V 1 

dx 

jl - (s?Ip 0 % 2 )_ 


f 

- (g 2 / Po 2c o 2 ) 


where 


g 2 = t X 2 + ty 2 


( 21 ) 
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Equation (20) implies the existence of a function 0 such that 


\jl ~ (g 2 / P 0 2 c q 2 ) 


i 1 - (^/Po^o 2 ) 


( 22 ) 


It is now readily seen that the above development is completely equiv- 
alent to the Chaplygin theory. However, the usual physical interpre- 
tations of 0 and the partial derivatives of 0 and are not 
implied. The formal relations, however, remain unchanged. 

It is well-known that the function 0 satisfies the equation for 
the minimal surface and that the solution for 0 and + can he 
expressed as 

o 


t - i0 = 2p 0 c 0 h(fi) 


(23) 


dz = dx + i dy 


dh 


= — + ft dh 

n 


(24) 


where h(ft) is an arbitrary analytic function, and 


n = 


% - 1% 


fpQ 2c o 2 - e 2 + Po c c 


(25) 


To satisfy the initial conditions (17), it is found convenient to apply 
the theory for minimal surfaces developed by Weierstrass, Bjorling, 
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and Schwarz (see reference ll). The final solution may he given in the 
following differential form: 


dx ■■ He 


F’ ( Pl UP) + i 


( x v 2 ) + k 2 ( f ') £ X + (F') g - Mf 
{ 1 + (F') 2 1 + (p-)2 + 2Lni 


5 «x 


(l - k 2 ) + k 2 & + ( F ') 2 - + (F*) 2 + 7%f] 


1 + (F 1 ) 2 


jl + (F' ) 2 + ~ • - M x 2 ‘ 


dy ■= Re 


1+1 


F' 

, . ^ i 1 " <‘" > 2 - “i 2 


1+ < F ') 2 X+ ( F .)2 + 7^_±m x 2 

I (l ' k 2 ) 

p (pi )2 [l + (F* ) 2 - Mj^J Jl + (F* ) 2 + 7M X ? 

1 + (! ’> 2 


> dP 


/ 




d0 + 1 d)jr 
P x u 


1 + (F ( ) 2 + 7 " , 1 Mi 2 


1 


+ 1 


(l . k 2 ) + k 2 ( F ') 2 E + (F’) g -M 1 2 J[l+ (F ') 2 + y Mx g 

1 + (F') 2 r ^ y - 1 ~I 2 

[i+ (F-) 2 + 2_J: Ml 2 


> dp 


M2 6 ) 


where Re denotes the real part of a complex quantity. 


k = Pl U/p 0 c 0 . 


and 


F' = F' (p-jUP) 

and P is a complex running variable. It may he noted that this solu- 
tion involves only the function F^p^y) for the slope of the shock 
wave. N ‘ 
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Weierstrass 1 formula for minimal surfaces .- Weierstrass* formula 
for minimal surfaces (reference ll) reads as follows: 



where f-^ is an analytic function of the complex variable t. 

Now 

t | ^’"(t) dt = dh (l/it) (28) 


From equation (27) the following relation is immediately derived: 


dz=dx+idy=itdh (l/it) + dh (l/it)/it 


Since dw = Re itf;]/" (t) dt = Re 2i dh (ft) = d0/p o c o , the above rela- 
tion reduces to equation (24) if one sets ft = l/it. 


Shock conditions and Bjorling problem .- As far as differential 
equation (19) is concerned, the parametric function h(ft) in expres- 
sion (24) always generates an exact solution. There are still shock 
conditions to be satisfied. These conditions (17) raise the problem of 
finding the solution in the large of the Cauchy initial-value problem 
for minima l surfaces. Geometrically, it is now reduced to the problem 
of finding the minimal surface- when a piece of space curve and the 
tangent plane at every point bn the curve are prescribed. The curve is 

given in T(x,y,w) space with the following parametric representation 
in vector form 



F(piUy) 


P X U 


y, 0/p o c o 




(29) 
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where 


0/p o c o = ^ 


F(plUy) 


Pl u 


= J 1 far (piUy) + ^JJiy 


with 0 X and 0y determined -by equations (22) and (17). The 
direction of the normal of the tangent plane at a point 

[*(pi u y)/pi u , y> JL] is given by (0 x ,0y,-p o c o ). This geometrical 

problem was first suggested by Bj or ling in 1844 and solved by Schwarz 
in 1874. . 


Schwarz formula .- Let £ be a vector with components x, y, 

and w. Let rf(P) and f^( Q) be complex- valued vector functions of 
the complex variables P and Q, respectively. The Schwarz formula for 

the Bjorling problem with the given curve £ = £ Q (r) and the normal of 

->2 

the tangent plane at any point on the curve k(t) such that k =1 
and k x £ 0 =0, where is the tangent vector of the curve, is 

( reference ll) : 


where 


21 = ^(P) + jf(Q) 

dif(P) = d'f Q (P) + ijjf(P) x 
dt(Q) = dt 0 (Q) - i[k(q) x 



(30) 

a?o( p D 

(31) 

<(«>] 

(32) 


From the last two relations it follows immediately that both vectors 
are isotropic lines; that is. 
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As the minimal surface should real one has only to take Q as the 
complex conjugate of P and £, as the complex conjugate vector 
function of T|. 

Exact solution of simplified differential equation satisfying 
shock conditions.- Now since 


K = 


*y/po c o> -^x/p o c o> “V 1 - (g 2 / Po 2c o 2 ) 


and ty^F' + = p-jU on the shock wave, it follows that 


K X d| Q = 


1 - Pl Ut y /p 0 2 c 0 2 -F’ + PjU^/p, 


o 2c o 2 Pl U 


f 1 " (^/Pq^O 2 ) j 1 ~ (fl P 0 2 °0 2 ) 


Pn c 


O^O 


dy 


= Jt(y) dy 


(33) 


where and are assigned values according to equation (17). 

Let the components of this vector he denoted hy Jjtj(y), rt 2 (y)> Jt^yOj dy. 

If F(p x Uy) is an analytic function of y and possesses an ana- 
lytic extension in a domain in the d-plane, where cr = y + its, then 
equation ( 30 ) is the differential form 


2 dft = Re [d? 0 (p) + ijf(P)j dP 


(3*0 


Or, 


dx = Re 


r 

If 


■ {^) 


q 2 q 2 


1 - Pl u* y / 

\ 1 - (s 2 |po 2c o 2 ) 


dP 


(35) 
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dy = Re 


1+ 1 - F ' (Pl OT ) + 

- (g 2 / Po 2 °o 2 ) 


(36) 


d0 = Re 


V(Pl°P) - *x 
^ - (s 2 /^ o 2 °o 2 ) 


+ ipjU 


dP 


(37) 


From the definition of h(fl), by equation (23), and relation (28), it 
is evident that the last equation in equations ( 27 ) can also he written 
as 


d 0 + i dijr 

Po c o 


itf*"(t) dt 


(38) 


Because of the fact that tj(P) in equation (31) is an isotropic line 
when ? = Constant, or (’l 1 ') 2 (P) + + = one can 

find an analytic relation which transforms t into P so that 


4 l‘(P) dP = i ^"'(t) dt ( 39 ) 

T) 2 '(P) dp = - ... — ^"’(t) dt (4o) 

and 

(41) 


t| 3 * (P) dP = itf 1 rt ' (t) dt 
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Therefore equation (37) can also be written as the following: 


d0 + i di}r = < 


[v , (pi up ) - *M - (e 2 /po 2c o 2 ) + iPlU 


^ dP 


(42) 


The final solution (26) is now obtained by substituting shock condi- 
tions (17) into equations (35), ( 38 ), .and (42). 


Massachusetts Institute of Technology 
Cambridge, Mass., July 6, 1949 
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